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Abstract 

Since  Pakistan  has  varying  climates  and  terrains ,  tlu  Pakistan  Arm;/  ro¬ 
tates  its  units  between  peacetime  locations  so  that  no  unit  endures  inequitable 
hardship  or  enjoys  unfair  advantage.  Army  policy  specifies  strict  constraints 
on  unit  rotations,  such  as  the  length  of  a  unit's  stay  in  any  location,  the 
number  of  units  moving  at  any  time,  and  the  allowable  replace  meats  for  any 
moving  unit.  Scheduling  rotations  manually  in  accordance  with  these  rules, 
as  is  currently  practiced,  is  extremely  difficult  and  lime-consuming.  This  pa¬ 
per  presents  an  integer  programming  model  that  finds  feasible,  minimum-cost 
schedules  for  the  Pakistan  Army's  desired  planning  horizons.  The  model  also 
ensures  that  the  units  arc  positioned  at  the  end  of  the  planning  horizon  so 
that  feasible  schedules  exist  for  future  planners.  The  model  is  implemented 
with  commercially  available  optimization  software.  Schedules  are  obtained 
for  realistic  test  problems  in  less  than  an  hour  on  a  personal  computer. 


The  Pakistan  Army  peacetime  rotation  problem,  as  described  in  the  ab¬ 
stract,  was  chosen  by  Colonel  Haig  as  a  research  topic  for  a  masters  thesis, 
Baig  [1992],  directed  by  his  co-authors.  lie  received  guidance  on  problem  def¬ 
inition  from  Brigadier  General  Maroon  Bashir  Khan,  who  was  then  serving  as 
the  Army  Attache  at  the  Pakistan  Embassy  in  Washington.  The  model  pre¬ 
sented  in  this  paper  captures  all  aspects  of  the  peacetime  rot  ation  problem 
as  described  by  the  general. 

The  peacetime  rotation  problem  is  a.  unique  scheduling/timelabling 
problem.  It  shares  the  diversity  and  large  number  of  constraints  typically 
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found  in  the  construction  of  sport  schedules.  (E.g.,  see  Andrue  and  Coromi- 
nas  [1989]  for  the  1992  Summer  Olympic  Games,  Shell  [1985]  for  the  Na¬ 
tional  Basketball  Association,  Cain  [1977]  and  Shell  [198.5]  for  Major  League 
Baseball,  and  Ferland  and  Fleurent  [1991]  for  the  National  Hockey  League). 
However,  the  peacetime  rotation  problem  has  relatively  few  events  compared 
to  the  hundreds  or  thousands  of  games  that  make  up  a  typical  sports  sched¬ 
ule.  An  exact  procedure  is  therefore  undertaken  in  this  paper  instead  of 
the  heuristic  procedures  used  successfully  for  large  scale  timetabling  prob¬ 
lems  Carter  [1986],  The  following  sections  present:  1)  policies  defining  the' 
problem,  2)  the  integer  programming  formulation,  3)  model  refinements  for 
increased  tractability,  1)  computational  experience,  and  5)  conclusions. 

1  Pakistan  Army  Peacetime  Rotation  Policies 

The  Pakistan  Army  classifies  military  locations  into  Peace  Areas  (PAs).  Semi- 
Hard  Areas  (SHAs)  and  Hard  Areas  (HAs).  This  classification  accounts 
for  Pakistan’s  diverse  terrain,  ranging  from  desert  to  lofty  mountains,  and 
temperatures,  ranging  from  below  freezing  to  above  10  degrees  Celsius.  The 
geographic  classifications  also  account  for  available  facilities  at  the  locations 
and  the  proximity  to  major  metropolitan  areas.  To  ensure  that  personnel 
servo  equally  in  all  three*  areas,  unit  personnel  are  rotated  between  locations 
during  peacetime  on  a  regular  basis. 

Military  units  in  the  Pakistan  Army  are  classified  in  three  ways:  according 
to  their  operational  role  (strike  or  defensive),  according  to  their  functional 
role  (Armor,  Artillery,  Engineers,  Infantry,  Signals,  Supply,  etc),  and  ac¬ 
cording  to  whether  their  equipment,  is  supplied  by  Eastern  or  Western  bloc 
nations.  (See  Cohen  [1981]  for  the  roles  of  the  East  and  West  in  equipping  the 
Pakistan  Army.)  Separating  the  units  into  distinct  categories  based  on  these 
three  attributes  simplifies  subsequent  mathematical  analysis  and  is  useful  for 
elucidating  rotation  policy. 
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The  General  Headquarters  of  the  Pakistan  Army  has  established  rotation 
policies,  whose  salient  features  are  as  follows. 

1.  A  unit  can  move  only  if  replaced  by  a  unit  of  the  same  category.  >.r., 
having  the  same  operational  and  functional  roles,  and  operating  the 
same  type  of  equipment.  This  results  in  a  natural  division  of  t  he  overall 
problem  into  separate  rotation  problems  for  each  category.  The  largest 
category,  currently,  has  87  units. 

2.  If  a  unit  at  location  A  moves  to  location  13,  then  a  unit  of  the  same 
category  at  location  J3  must  move  to  location  A  in  the  same  year. 
This  policy,  referred  to  as  mutual  replacement ,  simplifies  transfer  of 
operational  and  administrative  responsibilities  at  both  locations. 

3.  Equipment  does  not  move,  only  personnel  do. 

4.  No  more  than  one  unit  can  move  from  the  same  brigade  in  t  he  same 
year.  A  brigade  is  composed  of  three  units  of  the  same  category. 

5.  Some  units  (such  as  Engineers,  Signals,  Reconnaissance  and  Support 
Battalions)  do  not  belong  to  a  brigade.  They  fall  under  direct  control 
of  a  Division.  There  can  be  no  more  than  one  of  these  units  moving 
from  the  same  location  in  the  same  year. 

6.  Each  unit’s  tenure  requirement  varies  by  location  classification  as  fol¬ 
lows: 

•  Peace  Area  (PA)  5-7  years, 

•  Semi-Hard  Area  ( S 1 1 A )  2-4  years, 

•  Hard  Area  (HA)  1-3  years. 

7.  An  individual  unit  must  rotate  according  to  the  cycle  of  locations: 

PA  — »  S11A  — »  PA  — ►  11 A  — »  PA,  as  shown  in  Figure  1. 
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Figure  1:  Each  unit  must  rotate  through  Peace  Areas,  Semi-Hard  Areas  and 
Hard  Areas  in  the  indicated  order.  Alternate  Peace  Anas  in  the  cycle  must 
be  different.  The  length  of  slay  in  each  location  must  fall  with  in  a  prescribed 
range. 

8.  There  is  no  restriction  placed  on  which  SHA  or  HA  locations  a  unit 
visits  on  its  cycle.  However,  a  unit  must  not  return  to  its  previous  PA 
location. 

Currently,  peacetime  rotation  schedules  are  developed  manually  by  plan¬ 
ners  at  the  General  Headquarters  of  the  Pakistan  Army.  A  five-year  schedule 
is  developed  on  a  yearly  basis.  Due  to  the  large  problem  size  and  complex 
policy  structure,  the  current  method  suffers  from  the  following  drawbacks: 

•  The  units  may  be  positioned  at  the  end  of  the  planning  horizon  in  a  way 
that  precludes  future  schedules  from  satisfying  the  policy  constraints. 

•  The  schedule  requires  hundreds  of  man-hours  to  develop. 

•  It  is  difficult  to  evaluate  proposed  policy  changes. 

•  The  schedule  may  not  be  developed  impartially. 
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•  The  schedule  may  incur  excess  cost  by  transferring  units  more  often  or 
over  greater  distances  than  necessary. 

These  limitations  motivated  the  development  of  an  integer  programming 
model  to  assist  with  rotation  scheduling. 

2  Peacetime  Rotation  Model 

Our  integer  programming  model  for  the  Pakistan  Army’s  peacetime  rotation 
problem  minimizes  the  total  cost  associated  with  all  scheduled  moves,  while 
ensuring  that  all  policy  constraints  are  satisfied.  The  model  is  valid  for 

time  horizons  up  to  Ifi  years,  which  satisfies  the  Pakistan  Army’s  planning 

requirement  for  five-year  rotation  schedules. 

•  Indices: 

*,»'  units, 

/,  /'  locations, 

t,t'  years. 

•  Given  Data: 

iniili  the  initial  location  for  unit 

mini  minimum  stay  allowed  at  location  /, 

maxi  maximum  stay  allowed  at  location  l, 

mcosUw  movement  cost  for  unit  i  from  l  to 

stay n  number  of  years  unit  i  has  been  at  location  /  at.  the  start,  of  the 

first  time  period, 

V  set  of  peace  area  locations, 

‘H  set  of  hard  area  locations, 


S  set  of  semi- hard  area  locations. 

•  Derived  Sets:  The  following  sets  are  used  to  enforce  rotation  policies. 

Their  derivation  is  discussed  in  section  T 

Um  set  of  all  units  eligible  to  move  from  location  /  to  /'  in  year  /, 

F,n  set  of  all  possible  locations  from  winch  unit  i  could  have  moved  if  • 

it  arrives  at  location  /  in  year  /, 

Tut  set  of  all  possible  locations  to  which  unit  i  can  move-  if  situated  at 
location  /  in  year  /. 

•  Decision  Variables: 

Xui’t  1  if  unit  i  moves  from  location  l  to  location  /'  in  year  /,  and  0 
otherwise. 

•  Formulation: 


minimize  £  J2  52  vu  ostiW:r,w, 

i  1  i'  i 

mari  —  atai/ii 

52  52  =  1  V'1  1  =  unll< 

(0 

52  •r*f'n =  52  r'ii't  'vh'i 

•etVu 

(2) 

52  52  *;«'<  <  i  v//. 

v  ief-v. 

(3) 

i  i+lflftTi 

52  52  -  52  52  *»«'«' 

t'=i  »'€/',„<  t'=I+min, 

('!) 

f— mini 

XI  T»H't  $  XI  XI 

(•r>) 

t'eT.i,  t'=l-mar, 
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ZEL  -  1  v' 

/gr/'en  i 

(«) 

EEE  ^  1  v? 

ler/'e.s  < 

(7) 

EEE  r'«'<  ^  1  V' 

i^Hi'er  < 

(«) 

EEE  x'“'1  -  1  v' 

its  i'er  i 

(!») 

y  y  -I’, i'ii  +  y  y  <  j  v?,  /  e  v 

l'£H  i  t'£S  I 

(10) 

yyyy  = e  e  l  z  r>»’> 

•  leHi'er  i  i  ttsi'er  t 

(ID 

Constraint.  (I)  ensures  each  unit’s  first  move  is  completed  within  tin*  min¬ 
imum  and  maximum  toimrr  requirements  at  its  initial  location.  Constraint 
(2)  enforces  the  mutual  replacement  policy.  Constraint  (3)  allows  no  more 
than  one  unit  to  move  from  the  same  location  in  the  same  year.  If  two 
or  more  brigades  have  the  same  geographical  location,  this  location  is  bro¬ 
ken  into  separate  brigade  locations.  Thus,  constraint  (3)  also  enforces  t  he 
restriction  that  no  more  than  one  unit  can  move  from  the  same  brigade. 

Constraints  (4)  and  (-•>)  control  the  tenure  requirements  for  all  moves  tak¬ 
ing  place  after  the  first  move.  Constraint  (4)  restricts  the  number  of  times 
a  unit  enters  location  /  up  to  year  I  to  he  less  than  or  equal  to  t  he  number 
of  times  it.  leaves  the  location  up  to  year  /  -f  macj.  These  constraints  ensure 
that  a  unit  leaves  a  location  in  year  /  only  if  it  arrived  at  that  location  in  an 
appropriate  earlier  year.  (They  can  alternatively  be  formulated  with  noneu- 
nnilative  inequalities,  but  the  cumulative  form  provides  better  comput  ational 
performance.)  Constraint  (5)  requires  a  unit  leaving  location  /  in  year  1  to 
have  arrived  there  t  —  maxi  to  /  —  mini  years  earlier.  Constraints  (1)  and 
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(5)  cannot  prevent,  a  unit  from  moving  away  from  I  lie  same  locat  ion  more 
Ilian  once  in  the  last  few  years.  Constraints  (6)  -  (ft)  eliminate  this  problem 
by  allowing  each  unit  at  most,  one  move  between  different  areas  during  the 
planning  horizon  (a  valid  restriction  for  planning  horizons  of  10  years  or  less). 

Constraint  (10)  prevents  a  unit  from  moving  to  the  same  PA  more  than 
once.  If  the  planning  horizon  is  less  than  seven  years,  this  constraint  is 
unnecessary. 

Constraint  (1 1)  states  the  total  number  of  moves  from  II  As  to  PAs  should 
be  equal  to  the  total  number  of  moves  from  SI  I  As  to  PAs.  This  constraint 
helps  position  units  appropriately  at  the  end  of  the  planning  horizon,  as 
discussed  further  in  section  3. 

3  Model  Refinement 

Even  though  we  can  treat  the  Pakistan  Army's  rotation  scheduling  problem 
with  separate  models  for  each  operational/functional  area,  the  independent 
models  can  still  be  quite  large.  The  largest  operatiounl/fum  tional  area  lias 
87  units  spread  over  30  locations.  Straightforward  application  of  the  pre¬ 
ceding  model  for  these  units  over  an  eight-year  horizon  would  require  over 
600,000  binary  variables.  Fortunately,  characteristics  of  the  policies  can  be 
exploited  to  identify  many  impossible  unit  movements  and  eliminate  the  cor¬ 
responding  variables.  The  transition  eligibility  parameter  OKm't  is  defined 
for  this  purpose.  It  has  value  1  if  and  only  if  unit  i  is  eligible  to  move  from 
location  /  to  l'  in  year  t.  The  idea  is  to  make  this  parameter  zero  as  often  as 
possible  without  sacrificing  model  fidelity  or  optimality.  The  sols 
Fut  defined  in  section  2  are  readily  generated  from  this  parameter.  The  logic 
for  deriving  OK  mu  is  as  follows. 

1.  Rotation  policies  and  tenure  limits  severely  rest  rict  the  set  of  allowable 
moves.  For  example,  a  unit  that  has  moved  from  a  HA  and  has  been 
at  a  PA  for  three  years  can  not  possibly  move  for  two  more  years  and 
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then  only  to  a  S11A.  Subsequent  moves  are  similarly  restricted.  These 
situations  are  detailed  below  in  four  exhaustive  cases.  Lei  OP  A,  denote 
the  current  or  most  recent  PA  location  for  unit  i  in  the  initial  year, 

Stay  =  m\n{stayn,  mini)  for  the  unit’s  initial  locat  ion  /,  and  let  “\” 
denote  set  exclusion. 

a)  If  unit  i  is  currently  located  at  a  PA  and  previously  served  in  a  S 1 1 A , 

I  if  leV,  /'  6  H,  f,  <  /  +  Stay  <  8, 

l  if  leH ,  l'  e{V\  OP  A,},  7  <  M  Stay  <11. 

1  if  l  €{V\  0/M,},  l'  eS.  12  <  /  +  Stay  <  18, 

0  ol In  riris( . 

b)  If  unit  i  is  currently  located  at  a  PA  and  previously  served  in  a  IIA, 

'  1  if  le  v .  /'  e  5,  f>  <  /  +  stay  <  8. 

0K  t  =  i  1  if  i  e  S,  V  e{V\  OVA,}.  8  <  /  +  Stay  <  12, 

Xiirt  "|  i  if  i  e  {v  \  or  a,},  r  en,  1.3  <t  +  stay  <  19, 

(  0  ofhcru'is c. 

c)  If  unit  i  is  current  ly  located  at  a  SI1A, 

1  if  leS,  l'e{V\OrAi},  3  <l  +  Slay<r>, 

1  if  l  e  {V\  O/M,},  I'eH,  8  <  /  +  Stay  <  12, 

1  if  l  eH,  l'  €  V,  9  <  /  +  Stay  <  15, 

0  otherwise. 
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d)  If  unit  i  is  currently  located  at  a  II A . 


i  tf  leh,  v  e{V\  or  A,}.  2  <  *  +  stay  <  1, 

1  if  i  e  {v  \  or  a,},  v  e  s,  7  <  /  4  stay  <11, 

I  if  l  e  S,  V  e  V,  0  <  t  4  stay  <  15. 

0  otherwise. 

For  a  detailed  example,  suppose  unit  i  is  currently  located  at  a  PA  and 
its  old  location  is  a  II A  (case  b).  This  unit  can  move  to  a  SI  1 A  after 
completing  5  to  7  years  of  stay  at  the  PA,  so  it  is  eligible  to  move  only 
when  t  +  Stay  is  6,  7  or  8  years.  The  unit’s  next  move,  to  any  PA 
except  Or  A,  must  take  place  2  to  1  years  later,  /.  c.,  in  one  of  the  years 
8  through  12. 

2.  The  mutual  rotation  policy  also  helps  eliminate  many  variables.  Con¬ 
sider  a  unit  at  location  /  that  is  eligible  to  move  to  V  in  year  /  after 
executing  step  1.  That  move  can  be  scheduled  only  if  another  unit  is  el¬ 
igible  to  move  from  location  /'  to  location  /  in  the  same  year.  Therefore, 
any  OKavt  that  was  1  after  step  1  is  changed  to  0  unless: 

T.  OKi'iut  >  1. 

3.  We  can  extend  the  idea  of  step  2  to  subsequent  moves.  For  unit  i  to 
be  eligible  for  a  move  to  location  /'  in  year  t ,  there  must  be  another 
unit  i'  eligible  to  replace  unit  i  at  location  V  between  /  4  »nm/«  and 
t  4  mazy.  Therefore,  any  OK  mu  that  remains  I  after  the  first  two 
steps  is  changed  to  0  unless: 

l+mniji 

mi  n  oKmv «•  >  i  • 

t'/i  l  l'=1+min  (r 
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Figure  2:  A  feasible  six-unit  rotation  pattern  that  can  be  extended  indefinitely. 
Nodes  in  the  figure  represent  units  undergoing  moves  in  the  indicated  year. 
Edges  represent  mutual  replacement. 


3.1  Conditions  for  Future  Feasibility 

The  Pakistan  Army  needs  to  ensure  that  units  are  positioned  at  the  end 
of  the  planning  horizon  so  that  feasible  schedules  exist  for  future  planners. 
Sufficient  conditions  are  developed  for  this  purpose.  These  conditions  are 
explained  with  the  help  of  Figure  2. 

Figure  2  shows  an  indefinitely  repeatable  rotation  schedule  that  satisfies 
all  restrictions  of  the  rotation  policy  outlined  in  section  1.  This  figure  con¬ 
tains  4  PAs,  1  HA,  and  1  SHA  locations  where  each  location  has  one  unit. 
The  circles  contain  unit  identifiers  and  the  arcs  between  the  circles  indicate 
units  exchanging  locations.  For  the  units  initially  located  at  PAs  (1,  2.  o, 
and  6),  their  last  area  is  shown  with  H  (for  HA)  or  S  (for  SHA). 
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It  is  possible  to  add  units  in  multiples  of  6  (4  PA  units,  1  I1A  unit,  1  SI  I A 
unit)  up  to  the  maximum  of  three  at  each  location.  Hath  G  additional  units 
(1  at  each  existing  location)  could  be  feasibly  added  to  the  above  schedule 
by  allowing  moves  in  similar  3  year  increments  starting  at  year  1  for  the  first 
6  units  and  at  year  2  for  the  second  6.  New  locations  with  units  following 
the  same  pattern  could  also  be  added. 

The  positioning  of  units  in  f  igure  2  obeys  the  following  conditions  at  all 
times. 

Condition  1  The  total  number  of  units  in  PAs  is  twice  the  number  of  units 
in  II As  and  SII As. 

Condition  2  The  number  of  units  in  HAs  and  SIIAs  is  equal. 

Condition  3  Half  the  units  in  PAs  previously  served  in  SIIAs  and  half  pre¬ 
viously  served  in  llAs. 

The  three  conditions  above  are  not  necessary  to  guarantee  t  he  existence 
of  a  feasible  solution  in  the  future.  However,  as  proven  in  the  appendix,  if 
the  conditions  are  satisfied  for  the  six  years  preceding  the  current  horizon 
and  are  enforced  throughout  the  horizon,  then  future  feasibility  is  guaranteed. 
Unfortunately,  some  operational/functional  areas  have  not  always  conformed 
to  the  conditions.  Therefore,  the  test  problems  of  section  4  are  solved  wit  h 
constraint  (11)  instead  of  the  more  restrictive  form  of  t  he  constraint 

EEE r»«'« =  53 13 13 x'11''  v  ( 1  2) 

i  tewt'er  «  lesrer 

which  would  guarantee  feasible  rotations  in  the  future  with  appropriate  start¬ 
ing  conditions.  Though  not  guaranteed,  constraint  (11)  empirically  yielded 
rotation  schedules  with  ending  conditions  that  allowed  future  feasible  rota¬ 
tions. 
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4  Computational  Experience 

The  Pakistan  Army  peacetime  rotation  model  is  implemented  in  the  (Seneval 
Algebraic  Modeling  System,  GAMS,  [1992]  and  solved  using  XA  [1987]  and 
OSL  [1991].  All  computational  results  are  obtained  using  a  486/33  personal 
computer  with  16  megabytes  of  RAM.  The  goal  of  our  work  was  to  develop 
an  implementable  scheduling  model.  We  believe  our  choice  of  commercially 
available  software  and  a  personal  computer  represents  the  best  chance  of 
implementation  for  the  following  reasons: 

1.  The  Pakistan  Army  can  implement  the  model  for  a  reasonable  cost. 

2.  The  software  is  stable,  well  documented,  and  the  user  can  benefit  from 
future  software  improvements. 

3.  Software  maintenance  cost  is  low. 

4.  The  software  is  portable  to  new  platforms  and  operating  systems. 

•r>.  Algebraic  modeling  languages  such  as  CAMS  allow  easy  modification 
and  addition  of  constraints. 

We  develop  eight  test  problems  as  described  in  Table  1.  Problem  names 
are  vague  and  data  is  hypothetical  for  security  reasons,  but  problem  sizes 
are  representative  of  actual  Pakistan  Army  situations.  The  computing  time 
required  for  GAMS  to  generate  test  problems  is  given  in  'fable  I.  A  large 
portion  of  the  GAMS  time  was  taken  deriving  the  parameter  Ol\,uu- 

All  the  test  problems  satisfy  Condition  3  of  Section  3  for  the  initial  time 
period.  All  but  INFANTRY  I  also  satisfy  Conditions  1  and  2.  The  test, 
problems  are  all  feasible,  experimentation  with  other  starting  conditions 
often  resulted  in  infeasibility. 

The  test  problems  of  'fable  1  were  solved  for  integer  solutions  with  a 
10%  optimality  tolerance  (f.c.,  termination  occurs  when  the  first  solution 
guaranteed  within  10%  of  optimal  is  obtained),  using  both  the  XA  and  OSL 
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Problem  Type 


Problem  Size 


(.JAMS  Generation 
Time  (seconds) 


I  -  ■  — -  ■  -L 


INFANTRY  1 

6  years 

87  units 

30  locations 

2,562  constraints 

12,241  binary  variables 
74,571  nonzeros 

530 

INFANTRY  2 

6  years 

72  units 

24  locations 

1,599  constraints 

5,746  binary  variables 
35,037  nonzeros 

236 

ARTILLERY 

6  years 

54  units 

21  locations 

1,088  constraints 

3, 140  binary  variables 
18,530  non zeros 

154 

7  years 

54  units 

21  locations 

1,563  constraints 

5,263  binary  variables 
31,967  nonzeros 

200 

8  years 

54  units 

21  locations 

2, 135  constraints 

7,828  binary  variables 
49,393  nonzeros 

258 

ENGINEERS 

6  years 

36  units 

19  locations 

913  constraints 

2,193  binary  variables 
14,829  nonzeros 

89 

7  years 

36  units 

19  locations 

1,253  constraints 

3,421  binary  variables 
23,379  nonzeros 

116 

8  years 

36  units 

19  locations 

1,684  constraints 

4,949  binary  variables 
35,430  non zeros 

147 

Table  1:  Test  pivblnn  description  and  model  generation  time  on  a  j. 
personal  computer. 
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Problem  Type  Solution  Time 
Horizon  (Year)  seconds  (solver) 


Iterations 


3, 575  (XA) 
2,338  (OSL) 

37,919  (XA) 
7,002  (OSL) 

mmsikimm 

1,508  (XA) 

-152  (OSL) 

21,647  (XA) 
3,178  (OSL) 

95  (XA) 

173  (OSL) 

2,519  (XA) 
2,190  (OSL) 

t  =  7 

079  (XA) 

211  (OSL) 

14,909  (XA) 
1,998  (OSL) 

t  =  8 

oo  (XA) 

455  (OSL) 

oo  (XA) 
3,642  (OSL) 

105  (XA) 

105  (OSL) 

4,735  (XA) 
1,394  (OSL) 

t  =  7 

2,855  (XA) 

341  (OSL) 

54,709  (XA) 
3,381  (OSL) 

t  =  8 

8,950  (XA) 
1,027  (OSL) 

105,918  (XA) 
11,212  (OSL) 

Table  2:  Test  Problem  Solution  Times 


solvers.  All  tuning  parameters  wore  left  at  their  default  values.  Table  2 
demonstrates  these  results  and  highlights  OSL’s  superior  performance  in  all 
but  one  case.  XA  was  unable  to  solve  the  linear  programming  relaxation  of 
the  8-year  Artillery  problem  apparently  due  to  cycling. 

For  completeness  of  the  computational  testing,  all  problems  were  also 
solved  to  optimality  (i.e.,  with  the  optimality  tolerance  set  to  zero).  The 
integrality  gap  (difference  between  the  optimal  linear  and  integer  program¬ 
ming  objective  function  values)  is  zero  in  five  of  the  eight  problems  tested 
and  less  than  3%  in  the  other  cases.  For  the  five  examples  with  no  integrality 
gap,  solving  with  zero  tolerance  took  about  the  same  time  as  solving  with  a 
10%  tolerance;  but  the  other  three  examples  took  significantly  longer.  Since 
the  integrality  gap  was  so  small,  the  added  computational  time  for  the  zero 
tolerance  yielded  no  improvement  in  solution  quality.  Though  this  behavior 
cannot  be  guaranteed,  we  recommend  using  the  10%  optimality  tolerance 
(the  GAMS  default)  for  future  instances  of  this  problem. 

5  Conclusion 

Computational  experience  with  the  model  demonstrates  that  optimal  unit, 
rotations  can  be  developed  for  6-8  year  horizons  on  a  personal  computer  in 
less  than  an  hour,  using  the  GAMS  modeling  language  and  the  OSL  solver. 
The  Pakistan  Army  requires  live-year  rotation  schedules  for  planning  pur¬ 
poses,  so  the  6-8  year  schedules  we  have  obtained  are  more  than  adequate 
in  scope.  The  computation  times  are  considered  accept  able  and  represent,  a. 
significant,  improvement  over  the  hundreds  of  man-hours  current  ly  used  to 
solve  the  problem. 

This  work  offers  some  general  lessons  for  practical  application  of  opt  i¬ 
mization  modeling. 

•  Use  general  purpose  software  when  possible.  Compared  to  special- 
purpose  algorithms,  the  costs  of  development  and  long-term  mainte- 
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nance  are  much  lower.  Also,  general- purpose  solvers  can  much  more 
readily  adapt,  to  changes  in  the  problem. 

•  Using  an  algebraic  modeling  language  to  generate  the  model  facili¬ 
tates  rapid  assessment,  of  comput  at  ional  tract  abilit  y  and  allows  several 
solvers  to  be  tested  competitively  with  minimal  effort.  In  our  expe¬ 
rience,  integer  programming  applications  are  too  varied  for  any  one 
solver  to  always  outperform  the  others.  (E  g.,  though  XA  lost  to  OSL 
in  this  case,  it  has  been  our  preferred  solver  in  other  applications.) 

•  When  computational  effort  prior  to  optimization  is  devoted  to  the  elim¬ 
ination  of  unnecessary  variables,  this  refinement  can  often  make  large- 
scale  instances  of  real-world  problems  tractable. 

•  Mathematical  analysis  leading  to  additional  constraints  beyond  those 
specified  by  the  original  problem  statement,  such  as  the  constraints  on 
the  ending  conditions,  can  lead  to  better  solutions;  and,  in  the  rase  of 
integer  programming,  carefully  chosen  extra  const  raints  often  make  the 
model  easier  to  solve. 
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Appendix 

Theorem  If  Conditions  1,  2  and  3  of  Section  3  are  satisfied  for  six  years, 
they  are  sufficient  to  guarantee  a  future  feasible  rotation  schedule. 

PROOF: 

Condition  3  stipulates  that  half  the  units  in  PAs  must  have  previously 
served  in  II  As  and  half  in  SI  I  As.  This  status  is  maintained  on  a  yearly  basis 
provided  that  whenever  a  unit  moves  from  a  II A  to  a  PA,  another  unit  moves 
from  a  SHA  to  a  PA.  The  units  at  II As  and  SIIAs  can  therefore  be  separated 
into  pairings  which  satisfy  one  of  the  following  three  cases. 

•  CASE  1.  A  unit  has  been  at  a  II A  for  1  year  and  another  unit  lias  been 
at  a  SIIA  for  1  year. 

•  CASE  2.  A  unit  lias  been  at  a  II A  for  2  years  and  another  unit  lias 
been  at  a  SIIA  for  2  years. 

•  CASE  3.  A  unit  lias  been  at  a  HA  for  3  }'ears  and  another  unit  has 
been  at  SHA  for  3  years. 

It  is  shown  for  each  case  that  it  is  possible  for  the  units  at  the  II A  and 
SHA  to  feasibly  rotate  in  the  same  year  with  a  unit  currently  at  a  PA  location. 
Any  six  unit3  satisfying  Conditions  l,  2  and  3  for  the  last  C  years  and  feasibly 
rotating  can  be  used.  Without  loss  of  generality,  the  six  units  are  numbered 
according  to  Figure  2  (i.e.,  unit  3  is  at  IIA).  For  clarity,  the  unit  identifiers  are 
bold  faced.  Also,  as  in  Figure  2,  year  0  is  considered  the  first  year  available 
to  change  a  unit’s  location. 

CASE  I  ■ 

Unit  3  has  stayed  I  year  at  a  IIA.  'Phis  implies  unit  3  replaced  unit  0,  I 
year  ago  and  unit  0  replaced  unit  1  either  2,  3  or  4  years  ago. 

Unit  4  has  stayed  I  year  at  a  SIIA.  This  implies  unit  4  replaced  unit  5, 
I  year  ago  and  unit  5  replaced  unit  2  either  3,  4  or  5  years  ago. 
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Condition  3  ensures  that  each  year  a  move  from  a  II A  to  a  PA  occurs, 
a  move  from  a  S 1 1 A  to  a  PA  also  oc  c  urs.  Therefore,  we  need  only  consider 
unit  6  (5)  replacing  unit  1  (2)  cither  3  or  4  years  ago. 

These  conditions  provide  the  following  rotation  eligibilities: 

Unit  4  is  eligible  to  move  in  years  1,  2  or  3  and  Unit  1  is  eligible  to  replace 
unit  4  in  years  2,  3  or  4  if  unit  1  was  replaced  3  years  ago,  or  in  years  1,  2 
or  3  if  unit  1  was  replaced  4  years  ago. 

Unit  3  is  eligible  to  move  in  years  0,  1  or  2  and  Unit  2  is  eligible  to  replace 
unit  3  in  years  2,  3  or  4  if  unit  2  was  replaced  3  years  ago,  or  in  years  1,  2 
or  3  if  unit  2  was  replaced  1  years  ago. 

Therefore,  it  is  feasible  for  both  unit  3  and  unit  4  to  rotate  in  year  2. 
CASE  2 

Unit  3  has  stayed  2  years  at  a  11  A.  This  implies  unit  3  replaced  unit  6,  2 
years  ago  and  unit  6  replaced  unit  1  either  3,  4  or  5  years  ago. 

Unit  4  has  stayed  2  years  at  a  SIIA.  This  implies  unit  4  replaced  unit  5. 
2  years  ago  and  unit  5  replaced  unit  2  either  4,  5  or  G  years  ago. 

Condition  3  ensures  that  each  year  a  move  from  a  II A  to  a  PA  occurs, 
a  move  from  a  SIIA  to  a  PA  also  occurs.  Therefore,  we  need  only  consider 
unit  6  (5)  replacing  unit  1  (2)  either  4  or  5  years  ago. 

These  conditions  provide  the  following  rotation  eligibilities: 

i 

Unit  4  is  eligible  to  move  in  years  0,  I  or  2  and  Unit  1  is  eligible  to  replace 
unit  4  in  years  1,  2  or  3  if  unit  1  was  replaced  4  years  ago,  or  in  years  0,  1 
or  2  if  unit  1  was  replaced  5  years  ago. 

Unit  3  is  eligible  to  move  in  years  -1,  0  or  1  and  Unit  2  is  eligible  to  replace 
unit  3  in  years  1,  2  or  3  if  unit  2  was  replaced  4  years  ago,  or  in  years  0,  1 
or  2  if  unit  2  was  replaced  5  years  ago. 

Therefore,  it  is  feasible  for  both  unit  3  and  unit  4  to  rotate  in  year  I. 
CASE  3 

Unit  3  has  stayed  3  years  at  a  II A.  This  implies  unit  3  replaced  unit  6,  3 
years  ago  and  unit  0  replaced  unit  1  either  4,  5  or  G  years  ago. 

Unit  4  has  stayed  3  year  at  a  SIIA.  This  implies  unit  4  replaced  unit  5, 


3  years  ago  arid  unit  5  replaced  unit  2  either  5,  6  or  7  years  ago. 

Condition  3  ensures  that  each  year  a  move  from  a  II A  to  a  PA  occurs, 
a  move  from  a  SI  1 A  to  a  PA  also  occurs.  Therefore,  wo  need  only  consider 
unit  6  (5)  replacing  unit  1  (2)  either  5  or  G  years  ago. 

These  conditions  provide  the  following  rotation  eligibilities: 

Unit  4  is  eligible  to  move  in  years  -1,  0  or  1  and  Unit  1  is  eligible  to  replace 
unit  4  in  years  0,  1  or  2  if  unit  1  was  replaced  5  years  ago,  or  in  years  -1,0 
or  1  if  unit  1  was  replaced  G  years  ago. 

Unit  3  is  eligible  to  move  in  years  -2,-1  or  0  and  Unit  2  is  eligible  to  replace 
unit  3  in  years  0,  1  or  2  if  unit  2  was  replaced  5  years  ago,  or  in  years  -1,0 
or  1  if  unit  2  was  replaced  G  years  ago. 

Therefore,  it  is  feasible  for  both  unit  3  and  unit  4  to  rotate  in  year  0.  □ 
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